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Abstract 

A simple and natural Gao algorithm for decoding algebraic codes is described. Its relation to the 
Welch-Berlekamp and Euclidean algorithms is given. 
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ir~> . I. Introduction 

O 

In the recent article Gao [1] described a simple and natural algorithm for decoding algebraic 
codes in the class of algorithms decoding up to the designed error correcting capability. The 
asymptotic complexity of this algorithm coincides with the complexity of the best Reed-Solomon 

ft- 

decoding algorithms, and the description is the simplest for known algorithms descriptions. In 
this paper the Gao algorithm's relation to the Welch-Berlekamp [2] and Euclidean algorithms 
[3], [4] is given. 

II. Definitions and notations 

Let us define the (n, k, d) Reed-Solomon (RS) code over GF(g) with length n = q — 1, number 
of information symbols k, designed distance d = n — k + 1, q is prime power. 
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The RS code generator polynomial is 

b+d-2 

g(x)= J] (x~a l ), 

i=b 

where a is a primitive element of GF(g), b is any natural number. 
The received vector is represented as a polynomial 

n— 1 n—1 n—1 

R(x) = nx l = C(x) + E(x) = CiX* + £ e l x\ 

i=0 i=0 i=0 

where C(x) is the codeword, E(x) is the error vector. 

The ith error in the error vector E(x) has a locator Zi e {a , a 1 , a 2 , • • • , a"" 1 } and an error 
value Fi e GF(g)\0. 

The error locator polynomial is 

W(x) = f[(x-Zi), 
i=i 

where t < (d — l)/2, Zj is the locator of the error position in the error vector 

There are several coding methods for RS codes. In this paper two methods are used: the 
remainder coding in the time domain for systematic coding and the spectral coding in the 
frequency domain for non-systematic coding. Note that the decoding algorithm does not depend 
on the coding method. 

In Sections 4 and 5 remainder coding is applied. The codeword consists of two parts: the 
message part and the parity part. The message error locator polynomial is 

W m (x) = H(x-Z i ), 
i=i 

where t m < (d — l)/2, Z^ e a? , j e [d — 1, n — 1], is the locator of the error position in the 
message part of error vector E(x). 

In Sections 3 and 6 spectral coding is applied. The message polynomial of RS code is 

k-l 

M(x) = rriiX 1 . 

i=0 

The component q of codeword C(x) is computed as 

c i = M{a i ), ie[0,n-l]. 
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III. Gao algorithm 

We describe here the Gao algorithm [1] for the case of the classical RS codes only (n — q — 1). 

1) Interpolation. 

Construct an interpolating polynomial T(x) such that 

T{a l )= ri , ie[0,ra-l], 

where degT(x) < n. 

2) Partial GCD. 
Solve a congruence 

' W(x)T(x) = P(x) modx" - 1 
< degP(:r) < 2±* 
maximize deg P(x) 

by applying the extended Euclidean algorithm to x n — 1 and T(x), and we obtain unique 
pair of polynomials P(x) and W(x). 

3) Division. 

The message polynomial is 

MM = 

W(x) 

The asymptotic complexity of this algorithm 0(n(logn) 2 ) coincides with the complexity of 
good classical RS decoding algorithms [4], [5]. 

The first step of the Gao algorithm can be implemented by any fast algorithm for the discrete 
Fourier transform, for example, [6] . If the number of multiplications is to be minimized, then the 
best known algorithm of the discrete Fourier transform for small lengths (n < 511) is described 
in [7]. 

The best implementation of the second step is the Moenck algorithm [8], that is also reproduced 
in [9], [10]. Note that Moenck' s implementation of the second step completely coincides with 
the algorithm of the key equation solving by Sugiyama et al. [3]. 

For application of the Gao algorithm to other classes of algebraic codes, such as BCH, Goppa, 
or alternant codes, the additional re-encoding step is needed. 
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IV. The original Welch-Berlekamp algorithm 

The Welch-Berlekamp algorithm is the remainder decoding algorithm. We follow the original 
patent [2] and paper [11] with simple algorithm description. 

The Welch-Berlekamp algorithm consists of the following steps: 

1) Calculate the syndrome for cyclic codes 

d-2 

S(x) = ^ s,ix % = R(x) modg(x). 

2) Solve the key equation 

' Pj a?N(a?) = 8jW m {oP), j e[0,d- 2] 
degN(x) <degW m (x) < (d-l)/2, 



(1) 



where 

and obtain polynomials N(x) and W m (x). 

3) Determine the error locations in the message part of received vector from W m (x) and 
obtain the set {Zi} E a? , j E [d — 1, n — 1]. 

4) Determine the error values in the message part of received vector 

Y - HZ ) N{Zi) 
where W,^ n (x) is a formal derivative for W m (x), 

f{z) = z- b Y J Pl 



i=o a 



for Z E a j , j E [d - 1, n - 1]. 
5) Determine the parity part of codeword. 

V. Chambers interpretation 

In [12] Chambers described the second step interpretation of the Welch-Berlekamp algorithm. 
For RS codes pj ^ 0, j E [0, d — 2], and from the key equation (1) we have 

N{oP) = ^W m {pP\ je[0,d-2] 
degiV(x) < degW m (x) < (d - l)/2. 
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Construct an interpolating polynomial L(x) such that 



L(a j ) 




je[0,d-2], 



where degL(:r) < d — 1. 



From 



N(aj) = L(aP)W m (aP), je[0,d-2], 



we have 



' N(x) = L(x)W m (x) mod g(x) 
degN{x) < degW m {x) < (d - l)/2 
minimize deg W 7 ™ (x) . 



Solve this congruence by applying the extended Euclidean algorithm to g(x) and L(x), and we 
obtain polynomials N(x) and W m (x). 

VI. Another approach to the Welch-Berlekamp algorithm 

In [13] Gemmell and Sudan described a modification of the Welch-Berlekamp algorithm which 
is easy to understand. However, direct implementation of this modification is inefficient. 
If r, = Ci then ri = M(a i ). If r\ ^ a then W(a i ) = 0. Hence, 

= W{a i )M{a i ), i e [0, n - 1]. 

Let P(x) = W{x)M(x). Then the key equation is 



We solve the key equation and obtain polynomials P(x) and W(x). The message polynomial is 



The Gemmell-Sudan modification applies the idea of the original Welch-Berlekamp algorithm 
to the frequency domain. Clearly, the description of the Gemmell-Sudan modification is simpler 
than the description of the original Welch-Berlekamp algorithm. 



W(a i )r i = P(oi i ), ie [0,n-l]. 



(2) 
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VII. Improvement of the Gemmell-Sudan modification 

We apply the Chambers method to the key equation (2). 
Let us construct an interpolating polynomial T(x) such that 

T{a l )= ri , [0,n-l], 

where deg T(x) < n. This corresponds to the first step of the Gao algorithm. 
Further, from 

W(a i )T(a i ) = P(a% i e [0, n - 1], 

we have the congruence 

W(x)T(x) = P(x)modx n - 1. 

Solving this congruence corresponds to the second step of the Gao algorithm. 

It is clear that by combining methods from Sections 5 and 6, we get the Gao algorithm. 

VIII. Conclusion 

The paper does not contain any new decoding algorithms. However, the demonstrative con- 
nection between known decoding algorithms has important methodological significance in the 
coding theory. We see that the Gao algorithm belongs to the class of the frequency domain 
algorithms with simple description. Thus the paper demonstrates how carrying the algorithms 
from the time domain into the frequency domain results in simplification of their description. 
In author's opinion the Gao algorithm and its modifications are the simplest for codes with 
short lengths for any implementation. The correctness proofs of all considered algorithms are 
simple and can be found in corresponding references. The properties of the extended Euclidean 
algorithm are described in popular books on the coding theory [4], [14]. 
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